An analysis of the far-field response to external forcing of a suspension in Stokes fiow 

in a parallel-wall channel 

J. Blawzdziewicz 
Department of Mechanical Engineering, Yale University, P. O. Box 20-8286, New Haven, CT 06520 



00 

o 

O 
(N 

< 



o 



I 

O 

o 



> 

o 

^' 

o 

00 
O 



X 



E. Wajnryb 

IPPT, Swi§tokrzyska 21, Warsaw, Poland 

(Dated: August 13, 2008) 

The leading-order far-field scattered flow produced by a particle in a parallel-wall channel under 
creeping flow conditions has a form of the parabolic velocity fleld driven by a 2D dipolar pressure 
distribution. We show that in a system of hydrodynamically interacting particles, the pressure 
dipoles contribute to the macroscopic suspension flow in a similar way as the induced electric dipoles 
contribute to the electrostatic displacement field. Using this result we derive macroscopic equations 
governing suspension transport under the action of a lateral force, a lateral torque or a macroscopic 
pressure gradient in the channel. The matrix of linear transport coefficients in the constitutive 
relations linking the external forcing to the particle and fluid fluxes satisfles the Onsager reciprocal 
relation. The transport coefficients are evaluated for square and hexagonal periodic arrays of fixed 
and freely suspended particles, and a simple approximation in a Clausius-Mossotti form is proposed 
for the channel permeability coefficient. We also find explicit expressions for evaluating the periodic 
Green's functions for Stokes fiow between two parallel walls. 



I. INTRODUCTION 

There have been numerous studies of the effect of con- 
finement on the dynamics of rigid particles [1-8], de- 
formable drops [9-12], and macromolecules [13, 14] in 
creeping flow in a parallel-wall channel. The investi- 
gations revealed new confinement-induced phenomena 
such as migration of macromolecules away from the walls 
[13, 15-17], stability of strongly elongated drops in a con- 
fined shear flow [9, 10, 12], and cross-streamline migra- 
tion of spherical particles due to pair encounters in wall 
presence [18]. Confinement- related collective phenomena 
include spontaneous formation of string-like drop config- 
urations [9, 10], propagation of displacement waves in lin- 
ear trains of drops in Poiseuille flow [19, 20], instabilities 
of confined particle jets [21], and pattern formation and 
rearrangements of particle lattice in 2D regular particle 
arrays [20]. However, while the above studies revealed 
rich dynamics resulting from hydrodynamic interactions 
of particles with the channel walls, a comprehensive un- 
derstanding of mechanisms underlying the confinement 
effects is still lacking. 

Confinement-induced multiparticle collective phenom- 
ena often emerge as a result of hydrodynamic interactions 
associated with the far-field form of the flow produced by 
the particles moving in the channel [19-23]. Due to con- 
finement, this flow qualitatively differs from the far-field 
flow caused by particle motion in free space. The dif- 
ference stems from the strong fluid- volume conservation 
constraint associated with the wall presence, and from 
absorption of momentum by the walls. Owing to the 
momentum absorption, the velocity field decays too fast 
to produce a nonzero fluid flux through the boundary 
at infinity. Therefore, the fluid displaced by a moving 
particle creates a backflow pattern [4, 22-24] in order to 
ensure fluid incompressibility. In contrast, in unbounded 



systems the fluid in the whole space moves in the same 
direction as the particle. 

The far-fleld backflow produced by a particle in a chan- 
nel has a form of a parabolic Hele-Shaw flow driven by 
a 2D dipolar pressure distribution [4, 22-24]. The Hele- 
Shaw form of the flow far from the particle can be derived 
using an appropriate lubrication expansion [25]. The 
dipolar character of the velocity fleld around a spheri- 
cal particle results from the cylindrical symmetry of the 
problem. 

An immediate consequence of the backflow effect is 
the negative sign of the transverse component of the 
two-particle mutual hydrodynamic-mobility coefficient 
[3, 22]. The backflow also causes a large resistance of 
elongated particles in a narrow channel [3, 4, 7, 23]. Fur- 
thermore, collective action of the dipolar flow flelds pro- 
duced by individual particles gives rise to propagation of 
particle-displacement waves in linear particle arrays in 
Poiseuille flow [19], and it governs macroscopic deforma- 
tion and lattice rearrangements in regular particle arrays 
[20]. 

The dipolar far-fleld flow produced by the particles in a 
channel is crucial not only for understanding suspension 
dynamics on the particle scale, but also for describing the 
macroscopic fluid and particle transport. We show that 
particle contribution to the macroscopic volume flux in 
suspension flow through the channel can be determined 
from the amplitudes of the dipolar Hele-Shaw far-field 
flows produced by the particles. This behavior has a close 
physical analogy in electrostatics, where the induced elec- 
tric dipole moments of the molecules of a dielectric ma- 
terial contribute to the electric displacement field. We 
explore this analogy in our analysis of particle and fiuid 
transport through a channel. 

Our paper is organized as follows. Far-field dipolar 
particle response and particle polarizability are discussed 



in Sec. II, where we give a general outline of our theory. 
The dipole moment of an arbitrary induced-force distri- 
bution in a channel is evaluated in Sec. III. The polariz- 
ability cocfhcicnts for a spherical particle are determined 
in the friction and mobility formulations in Sec. IV. The 
results from the preceding sections are used to analyze 
macroscopic suspension flow in Sees. V and VI: Sec. V 
relates the macroscopic suspension velocity to the dipolar 
density per unit area of the channel, and Sec. VI provides 
macroscopic transport equations describing the particle 
and fluid transport. Our conclusions are drawn in Sec. 
VII. 

An important additional result of this study is the 
derivation of explicit Ewald-summation expressions for 
the periodic Green's functions for Stokes flow in the 
parallel-wall geometry. These formulations are used in 
our current numerical simulations to supplement the the- 
oretical analysis, and our expressions are also applicable 
in Stokesian-dynamics and boundary-integral algorithms 
for dispersion flows in parallel-wall channels. 



II. PARTICLE POLARIZABILITY 

In this section we focus on a general discussion of the 
far-field response of a spherical particle in a parallel- wall 
channel to an applied lateral force J- ^ torque I", and 
Poiseuille flow v'^'*^* driven by a constant lateral pres- 
sure gradient V||p™*. This response is represented in 
terms of particle polarizability coefficients, which are sub- 
sequently used to determine the macroscopic constitutive 
relations describing suspension transport in a channel. 

Our results for the particle response to external forc- 
ing can also be used to construct a simplified descrip- 
tion of particle dynamics in parallel-wall channels. In 
this simplified formulation (summarized in Sec. II C 2), 
the interparticle hydrodynamic interactions are incorpo- 
rated solely through the far-field flow, in analogy to the 
point-particle approximation for unbounded systems. As 
discussed in [20] , such a single-scattering approximation 
yields accurate results if the interparticle separation is 
sufficiently large. 



A. System definition 

We consider the dynamics of a spherical particle of 
radius a (or an array of such particles) in a parallel-wall 
channel of width _ff, under creeping-flow conditions. The 
walls are in the planes z ~ and z ~ H, and the particle 
center is at the axis z at a distance z = Z from the 
lower wall. The fluid velocity field satisfies the no-slip 
boundary conditions on the walls and the particle surface. 

For simplicity, the analysis presented in Sees. II-IV is 
for a single particle in the channel. However, our theo- 
retical formulation can be readily generalized to multi- 
particle systems. We use such a generalization in Sees. 



V and VI, where we consider macroscopic response of 
suspension to external forcing. 



B. Far-field scattered flow 

In the near-field regime p ^ H (where p = |p|, and 
p — xex + yGy is the lateral position with respect to 
the particle center) the flow fleld scattered by the par- 
ticle, v', has a complex 3D form that involves multiple 
image singularities [26]. However, in the far-field regime 
(yo 3> H) the scattered flow tends exponentially (on the 
lengthscale H) to a much simpler 2D Hele-Shaw flow of 
the form [23] 



^r'{r) = -^Tj-^z{H-z)V\\p'ip), 



(1) 



where rj is the fluid viscosity, p' is the perturbation pres- 
sure, V|| is the gradient operator with respect to the 
lateral coordinates p, and r = p + ze^- 

Away from the singularity at p = 0, the pressure field 
p' satisfies the 2D Laplace equation 



Vfp'(p) = 0, p ^ 0, 



(2) 



owing to the flow incompressibility. For a single particle 
moving under the action of an external force, torque, 
or Poiseuille flow, the pressure p' assumes the form of 
a 2D dipolar field, due to the cylindrical symmetry of 
the problem and the vectorial character of the forcing. 
Accordingly, we have 



p'(p) = -;1d.Vi,$-(p) = ^d4, 



where 



27r 



*o (p) = - Hp) 



27r 



is the solution of the 2D Poisson equation 

Vl%{p) = -27rS{p), 



(3) 



(4) 



(5) 



and D is the dipole moment of the perturbation pressure 

P'- . 

As discussed in the introduction, the Hele-Shaw flow 

(1) with the dipolar pressure distribution (3) involves the 

backflow effect. Namely, on the symmetry axis parallel to 

the dipole moment D, the fluid displaced by the particle 

is moving in the direction of D, whereas on the transverse 

axis it is moving in the opposite direction. 



C. Polarizability and mobility relations 

Since the dynamics of the system is governed by the 
linear Stokes equations, the dipole moment D is linear in 
the strength of the external forcing. 



T^'?" 



^H^ID = pP^-J^ + pF' ■ r + /xPPViip 



(6) 



(where the factor j^ri~^H^ is introduced to ensure the 
Lorentz symmetry for the matrix of transport coefficient, 
as explained in Sec. IV). The dipolar far- field perturba- 
tion pressure (3) is analogous to a 2D electrostatic dipo- 
lar potential. Relation (6) will thus be termed a polariz- 
ability relation, and the linear transport coefRcients fi^^ 
(A = p, t,r) will be referred to as the polarizability co- 
efficients (by analogy with an electrostatic problem of 
polarizable particles in an external electric field). 

Apart from producing the far-field response described 
by Eqs. (l)-(6), a particle in a channel also undergoes 
translational and rotational motion with the linear and 
angular velocities U and ft. This rigid-body particle mo- 
tion is characterized by the standard mobility relations 



U = fj,**.J^ + fj,*'.T-fj, 



tp, 



'\\P 



n^fi''-T + n"-r-n'p- v^p" 



(7a) 



(7b) 



where fi denotes hydrodynamic mobility coefficients. 

By a partial inversion of equations (6) and (7), the par- 
ticle response can also be expressed in a friction-relation 
form, where the velocities U and H are the independent 
variables, and the forces ^ and T are the dependent 
quantities (cf. Sec. IV). 

Equations (6) and (7) are related. Namely, as shown 
in Sec. IV, the matrix of mobility and polarizability co- 
efficients /i^^ satisfies the Lorentz symmetry 



/^ 



AB 



fJ^ 



BA-\ 



A,B = t,r,p. 



(8) 



where the dagger denotes the transpose of a tensor [27]. 
Equation (8) has significant consequences for suspension 
dynamics in parallel-wall channels, because it implies the 
corresponding Onsager reciprocal relation for the matrix 
of kinetic coefficients in the linear constitutive relation 
describing macroscopic suspension transport. 



1. Macroscopic suspension flow 

The polarizability relation (6) (and its multiparticle 
generalization) is essential for theoretical understanding 
of the macroscopic behavior of suspensions confined in 
a parallel-wall channel. This relation is also important 
for numerical evaluation of effective transport coefficients 
governing the suspension flux. 

In the following sections we demonstrate that the aver- 
age suspension velocity u can be expressed in terms of the 
average dipole moment of the particles D. Specifically, 
in Sec. V it is shown that 



u = Ko(-V||p + nsD) 



where 



Ko 



hr'H^ 



(9) 



(10) 



is the permeability of a particle-free channel, V||p is 
the macroscopic pressure gradient, and Ug is the particle 
number density per unit area of the channel wall. Equa- 
tion (9), supplemented with the polarizability relation 
(6) (averaged over the particle distribution) describes the 
macroscopic suspension flow (cf. Sees. V and VI). 

An intuition regarding physical interpretation of Eq. 
(9) can be gained by considering its electrostatic anal- 
ogy. A comparison of the hydrodynamic and electro- 
static problems indicates that the macroscopic pressure 
field is analogous to the electrostatic potential, and the 
pressure gradient corresponds to the electric field. The 
macroscopic suspension velocity u (which is divergence 
free due to the fluid incompressibility) corresponds to 
the electric-displacement fleld in the absence of external 
charges. Relation (9) is thus similar to the expression for 
the electrostatic displacement field in terms of the elec- 
tric field and the induced-dipole-moment density [28] . In 
Sec. V we further explore this analogy. 

A macroscopic theory based on the transport equa- 
tions (6), (7), and (9) is capable of describing complex 
phenomena that occur in suspension flows in a parallel- 
wall channels. In particular, we have shown in our recent 
study [20] that such a theory predicts a fingering insta- 
bility in evolving 2D particle arrays (which has been con- 
firmed by direct numerical simulations). 



2. Single- scattering approximation 

The amplitude D of the dipolar Hele Shaw flow scat- 
tered by a particle in a channel is a key quantity in a 
single-scattering approximation, where the one-particle 
polarizability and mobility relations (6) and (7) are com- 
bined with an assumption that the incident flow acting 
on a given particle is a superposition of the external flow 
and the far-field dipolar flows (1) produced by other par- 
ticles. 

The single-scattering approximation describes particle 
dynamics in dilute suspensions in parallel-wall channels 
under the strong-confinement condition H ^ 2a. Many 
fundamental collective phenomena in confined suspen- 
sion flows (e.g. instabilities of conflned particle jets [21], 
propagation of particle-displacement waves in ID arrays 
of drops [19], and pattern formation in 2D regular parti- 
cle arrays [20]) are driven by the far-field interparticle in- 
teractions associated with the dipolar scattered flow (1) 
and (3). Not only can essential qualitative features of 
such phenomena be captured using the single-scattering 
approximation, but it also gives accurate quantitative re- 
sults if the interparticle distances are sufficiently large 
[29]. 

Quantitative predictions using the single-scattering ap- 
proach require predetermination of the polarizability and 
mobility coefficients in Eqs. (6) and (7). We evaluate 
only the polarizability coefficients fi^^ (B = t,r,p), be- 
cause the mobility coefficients fi^^ (A = t, r , i? = t, r, p) 
for spherical particles between two parallel walls have al- 



ready been calculated [1~4, 23]. 



III. DIPOLE MOMENT OF INDUCED-FORCE 
DISTRIBUTION 

To evaluate the induced dipole moment of a particle in 
a channel, we need to solve the corresponding Stokes-flow 
problem. For this purpose, we apply the induced-force 
formulation and multipolar-expansion techniques. In this 
section we use the asymptotic Liron Mochon formula for 
the Green's function for Stokes flow in the parallel-wall 
geometry [24] , to relate the dipole moment D to the force 
distribution induced on the particle. We employ these re- 
sults in Sec. IV to calculate the polarizability coefficients, 
and in Sees. V and VI to determine constitutive relations 
describing macroscopic suspension flow. 



A. Induced-force formulation 

In our approach, the effect of a particle on the sur- 
rounding fluid is represented in terms of the induced- force 
distribution on the particle surface 



F(r)=a-25(ri-a)fi(ri), 



(11) 



where ri = r — Ze^ denotes the position with respect 
to the particle center, and ri = jrij. By definition of 
the induced force, the flow and pressure fields produced 
by the distribution (11) are identical to the velocity field 
v(r) and pressure p{r) in the particle presence [30-32]. 

For a particle in an external flow v°''*(r), the velocity 
and pressure fields can be represented by the boundary 
integrals 



v(r) =v'=^*(r)+ / T(r,r')-F(r')dr', 



p(r) = p'^^'ir) + / Q(r, r') • F(r') dr'. 



(12a) 



(12b) 



where p'^^^ (r) is the external pressure associated with the 
flow v'^'^*(r), and T and Q are the velocity and pressure 
Green's functions for a parallel-wall channel. 

The Green's functions for Stokes flow in the parallel- 
wall geometry were investigated using several differ- 
ent methods [1-4, 23, 24, 26]. In the Cartesian- 
representation approach proposed by our group [3, 4, 23], 
the velocity and pressure Green's functions T(r, r') and 
Q(r,r') are represented in terms of lateral Fourier inte- 
grals of simple matrix products. The explicit formulas 
are listed in Appendix A. 

For a particle moving with the translational and angu- 
lar velocities U and CI, the flow field (12a), evaluated at 
the particle surface S, equals the rigid body velocity of 
the particle 



With the above boundary condition, Eq. (12a) yields a 
boundary-integral equation for the induced forces. 

For a given induced- force distribution (11) the force 
and torque acting on the particle can be evaluated using 
expressions 



J^^ 



/ F(r) dr, T = / r X F(r) dr. 



(14) 



In the following section we determine the corresponding 
relation for the amplitude of the dipolar far-field flow 
produced by the particle. 



B. Dipole moment 

The dipolar strength D of a force distribution F can 
be obtained from Eqs. (12) using asymptotic far-field ex- 
pressions for the Green's function T and Q. (Such ex- 
pressions were first derived by Liron and Mochon [24] ; an 
alternative and much simpler derivation via a lubrication 
expansion is given in [25]). 

As discussed in [23], the asymptotic form of the 
flow and pressure Green's functions, THs(r,r') and 
QHs(r, r'), can be represented by the following formu- 
las 



THs(r,r') = -ir;-iz(if-z)VQHs(r,r'), 



QHs(r, 



-3tt-^H-^V%{p - p')z'{H 



(15a) 



(15b) 



where p = xe^ + ySy and p' = x'e^ + y'^y are the lateral 
position vectors, and $q"(p) is the point-source solution 
(4) of the 2D Poisson equation (5). We note that the 
asymptotic Hele-Shaw flow and pressure fields (15) sat- 
isfy the Stokes equations exactly. However, these fields 
have a different singularity than the original point-force 
singularity of the full Green's functions T and Q. 

According to Eqs. (15b) and (4), the far-field pressure 



/(P) 



QHs(r, 



•Fii 



(16) 



produced by a lateral point force F|| has a form of a 2D 
potential dipole 



/(P) 



1 

'2^ 



Do-V||$o(p-p': 



(17) 



v(r) 



rb 



(r)EEU + r2 X ri, res'. (13) 



with the dipole moment 

Do = 6H-^z'{H - z')F|| . (18) 

The streamlines of the corresponding velocity field 

vd(r) = -i,y-iz(i7-z)V||/(p) (19) 

also follow a 2D dipolar pattern. The quadratic depen- 
dence of the dipolar strength (18) on the position z' of the 
point where the force is applied follows from the Lorentz's 
symmetry of the Green's function (f5a) [23]. We note 



that the flow field produced by a transverse force (i.e. a 
force pointing in the z direction) is exponentially small 
in the far-field domain. 

Similar to the corresponding electrostatic problem, the 
dipole moment (18) represents not only the amplitude of 
the far-field pressure (16) but also the strength of the 
dipolar-pressure source, 



Vf/(p) = Do-V|,<5(p-p')- 



(20) 



In Sec. V the above expression will be used in our deriva- 
tion of the relation between the macroscopic dipolar- 
strength density per unit area of the channel and the 
average suspension fiow. 

The dipole moment of the far-field perturbation pres- 
sure (3) produced by a particle in a channel is evalu- 
ated by integrating (18) over the particle surface. Taking 



II 



l||F(r), where 



^.T^.T 



(21) 



denotes the projection operator onto the lateral direc- 
tions X and y, and F(r) is the induced-force distribution 
(11), we get 



D = 6iJ-3 /" z'(H _ z') l|| • F(r') dr'. 



(22) 



In the following section, Eq. (22) is used to determine the 
polarizability coefficients of a particle. 



IV. POLARIZABILITY COEFFICIENTS 

A. Multipolar expansion 

To determine the force distribution induced on the 
particle surface, the boundary-value problem (12a) and 
(13) is solved using the multipolar-expansion technique 
[3, 4, 23]. In our method, the induced-force distribution 
and the flow fleld in the system are expanded into the 
conjugate sets of basis functions introduced in [33]. In 
particular, we have the expansions 



F(r) 



E/('™'^)^w(ri) 



'(r) 



'(r) 






c(/TO(T)v+„^(ri), 



(23a) 



(23b) 



where the left-hand side of Eq. (23b) describes the exter- 
nal flow with respect to the rigid-body particle motion 
(13). 

In the above relations w^ and v^ are the conjugate 
spherical basis functions associated with the non-singular 
solutions of Stokes equations in spherical coordinates [3, 
4, 33], and f{lmu) and c{lrau) are the corresponding 
expansion coefficients [34]. The fiow v' scattered by the 
particle has a similar expansion in terms of singular basis 



functions V;~^^. Here 1 = 1,2,... and m = 0, ±1, . . . , ±1 
are the spherical harmonic orders of the basis functions, 
and the index a = 0, 1, 2 corresponds to the three types 
of Lamb's solutions for Stokes flow. 

By inserting expansions (23) into relation (12a) eval- 
uated at the particle surface one gets a linear algebraic 
equation of the form [3, 4] 



Y^ G{lma I l'm'a')f{l'm'a') = c{lma), 



(24) 



where the matrix elements G{lma \ I'm' a') are deflned in 
terms of the multipolar projections of the Green's func- 
tion T, 

G{lma I I'm' a') = (w+„,(ri) | T | w+„,^,(ri)). (25) 

In Eq. (25) we use the Dirac's bra-ket notation [35] 



(A IB) 



yA*(r).B(r)dr, 



(26) 



with the asterisk denoting the complex conjugate. 

Equation (25) is valid for a particle with the no-slip 
boundary conditions on its surface. For other boundary 
conditions (e.g., spherical viscous drops) the diagonal ele- 
ments of (25) need to be modified to include an appropri- 
ate single-particle scattering matrix [3, 4, 23]. It should 
also be noted that equation (24) is given here for a sin- 
gle particle, but general expressions for a multiparticle 
system are readily available [3, 4, 23]. 



B. Generalized friction matrix 

The response of the system to a given rigid-body par- 
ticle motion (13) and external parabolic flow 



-lr]-^z{H-z)Vp'' 



(27) 



:f 




- ^tt ^tr ^tp - 




u 


T 


= 


^rt ^rr ^rp 


• 


n 


_ ^i?3D _ 




^pt ^pr ^pp 




^-l^pcxt 



driven by a constant lateral pressure gradient Vp°''* can 
be characterized by the generalized resistance relation 



(28) 



As discussed below, the factor -p^H^ that multiplies the 
dipole moment D is needed to ensure the Lorentz sym- 
metry (35) of the resistance matrix C,. The coefficients in 
the bottom row of the matrix C, characterize the particle 
polarizability. 

The generalized resistance tensors C,^^ can be deter- 
mined by solving Eq. (24) to evaluate the induced-force 
amplitudes f{lma) in terms of the flow expansion coef- 
ficients c{lm(j). The solution can be represented by the 
matrix relation 



/(/too-) = ^ F{lma \ l'm'a')c{l'm'a'), 



(29) 



where F(lma \ I'm' a') denotes the elements of the matrix 
inverse to the Green's matrix G{lma \ I'm'a'). 

To obtain the resistance formula (28), Eq. (29) has to 
be supplemented with appropriate expressions that re- 
late the force ^, torque T , and dipole moment D to 
the expansion coefficients f{lma) of the induced force 
distribution (23a). We also need the relations for the ex- 
pansion coefficients c{lma) associated with a given rigid- 
body particle motion (13) and external flow (27). These 
relations can be expressed in terms of the transformation 
vectors X, discussed in Appendix B. Namely, we have 



and 



:F = ^X(t I lma)f{lma), 



(30a) 



cilma) = X(?mcr I t) • U -h X(/mcr I r) • n 

+ X{lma\p)-f]-^Vp''''\ (31) 



where the transformation vectors X satisfy the symmetry 
relation 



X(;mcr I A) =X*(A I Ima). 



(32) 



T = ^X(r I Ima) J {Ima), 



Ima 



1 tr3T~» _ \ ^ 



12 



H-'B 






X(p I lma)f{lm,a), 



Explicit expressions for the transformation vectors (32) 
(30b) aj-g listed in Appendix B. We note that only several 
matrix elements with small values of the indices /, m, a 
arc nonzero, according to Eqs. (Bl) and (B2). 

(30c) Combining Eqs. (29)-(31) yields the relation 



*/lB 



EE 



X{A I lma)F{lma \ l'm'a')X{l'm'a' \ B), A, B = t, r, p. 



(33) 



for the translational, rotational, and polarizability com- 
ponents of the generalized friction matrix. Taking into 
account relation (32) and the Lorcntz symmetry of the 
matrix F [3, 4], 



F{lma I I'm'a') = F*{l'm'a' \ Ima), 



(34) 



we obtain the corresponding symmetry of the generalized 
friction matrix 



ogous to (28) we have 



12' 



where 



u 




fi'^ 


/.^■^ 


VL 


— 


M^* 


M" 


-^H'^B 




/.P* 


/xP^ 






tp 






/•tt /-tr 



r 

_ypOxt 



(36) 



(37) 



-AB 



^BA\ 



(35) 



The symmetry relation (35) with A,B = t, r corre- 
sponds to the well-known Lorentz symmetry of the stan- 
dard resistance matrix [36]. For A = p or i? = p, how- 
ever, the symmetry relation (35) is new, and its impor- 
tance becomes clear in Sec. V, where we explain the re- 
lationship between the dipole moment induced on the 
particles and the macroscopic suspension flow. 



C. Mobility formulation 

The polarizability and mobility relations (6) and (7) 
can be obtained by partially inverting the generalized re- 
sistance relation (28). In the matrix representation anal- 



>*p ■ 


= 


■ Att Atr ■ 


-1 


■ ^tp ■ 

^rp 



[ /xP* flP' ] 



^2^pp = ^PP _ [ ^Pt ^pr ] 



tp 






/-tt 

/-rt 






(38) 



(39) 



Atp 

/•rp 



(40) 

Similar to the corresponding property of the friction ma- 
trix, the generalized mobility matrix /x satisfies the 
Lorentz symmetry (8), which is obtained by combining 
(35) with (37)-(40). Note that the negative sign is incor- 
porated into the pressure-gradient term on the right-hand 
side of Eq. (36) to obtain positive signs of the "pt" and 
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FIG. 1: Normalized polarizability coefficient (43), versus par- 
ticle position scaled by the vertical space available in the chan- 
nel, for channel width H/{2a) (as labeled). Coefflcient i^f^ 
characterizes the amplitude of the far-field flow produced by 
an immobile particle in an external pressure-driven flow. 
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FIG. 2: Normalized polarizability coefficient (45), versus par- 
ticle position scaled by the vertical space available in the chan- 
nel, for channel width H/{2a) (as labeled). Coefficient jl\'^ 
characterizes the amplitude of the far-field flow produced by 
a freely moving particle in an external pressure-driven flow. 



"pr" components of the matrix /x. We have also included 
factors ?7 and rp in Eqs. (38) and (40), to ensure that the 
polarizability and mobility coefficients have appropriate 
dimensionalities . 

The generalized mobility and friction tensors (28) and 
(36) are given here only for a single sphere. However, 
their multiparticle generalizations can easily be obtained. 



D. Tensor ial form 

For a single spherical particle, considered herein, the 
friction and mobility tensors i/)"^^ = C,"^^ , fi^^ are in- 
variant with respect to rotation around the axis z. Not- 
ing that T and Jl are pseudo- vectors, and that Vp'^'^* 
has only lateral components we find that 



ip 



AA 



-^^^1 



-V'f^e^e^, 



A = t, 



(41a) 



E. Numerical results 

In this section we present our numerical results for the 
generalized resistance and mobility matrices ^ and /x. 
Since we are concerned here with the effect of the far- 
field flow on the system dynamics, we focus on the polar- 
izability components that relate the external forcing to 
the induced dipole moment of the particle. 

We note that the translational and rotational compo- 
nents of the friction and mobility matrices were calcu- 
lated by Jones [2] and by our group [3, 4]. The effect 
of an external Poiseuille flow on particle motion has also 
been determined [1, 2, 23]. The polarizability coefficients 
■0P* and ■0'"^ (V' = Cj m) can be calculated from these ear- 
lier results by invoking the Lorentz symmetries (8) and 
(35). The coefficient ipPP, introduced here, has never 
been considered. 



1. Immobile particle 



^ 



AB 



^ 



AB 



A==p, B==t,p, 



i.AB 



(41b) 



,/,--=. ^--e||, yl = t,p, B = r, (41c) 

where L is the lateral unit tensor (21), and 

6|| Gx^y ^y^x V / 

is the lateral alternating tensor. The remaining compo- 
nents of the tensors xjj^^ are obtained using the recipro- 
cal relations (8) and (35). In Eqs. (41), the lateral trans- 
port coefficients are denoted by ijj , and the transverse 
coefficients by V'f'^- 



The dependence of the polarizability coefficient C^^ on 
the channel width and the transverse particle position Z 
is depicted in Fig. 1. The results are presented for the 
normalized coefficient 



cr 



pp _ 



12CP 



■K-qH^a^ 



(43) 



where the subscript 1 indicates that the quantity (43) is 
evaluated for an isolated particle in the channel. Accord- 
ing to equation (28), the coefficient (43) represents the 
dimensionless dipole moment 
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FIG. 3: Translational polarizability coefficient /i''', normal- 
ized by the corresponding result for the point force (47), ver- 
sus particle position scaled by the vertical space available in 
the channel, for channel width H/{2a) (as labeled). Coef- 
ficient /i^* characterizes the amplitude of the far-field flow 
produced by a particle moving under the action of a lateral 
force. 



FIG. 4: Translational polarizability coefficient /i'"^, normal- 
ized by the corresponding result for the point torque (48), 
versus particle position scaled by the vertical space available 
in the channel, for channel width H/{2a) (as labeled). Co- 
efficient /i'"^ characterizes the amplitude of the far-field fiow 
produced by a particle moving under the action of a lateral 
torque. 



produced by the imposed pressure gradient Vp"'^*. As 
discussed in Sec. VI A, the normalized polarizability co- 
efficient Cf appears as the leading-order term in the 
volume-fraction expansion of the permeability of a fixed 
bed of particles in the channel. The normalization has 
been chosen to simplify the resulting expression. 

The results in Fig. 1 indicate that the induced dipole 
moment (44) has a pronounced maximum for a particle in 
the midplane of the channel, especially for larger values 
of the channel width H . This behavior results from two 
factors: First, the flow acting on an immobile particle is 
the strongest in the center of the channel; and second, the 
induced force applied at point z' produces the strongest 
far-field flow for z' — H/2, according to Eq. (18). 



2. Suspended particle 



For a freely moving suspended particle, the numerical 
results are shown for three polarizability coefficients: jjP^, 
fi^*^ , and /iP'' . These coefficients characterize the induced 
dipole moment produced by the external flow, force and 
torque acting on the suspended particle. 

Figure 2 illustrates the dependence of the polarizabil- 
ity coefficient /i^P on the channel width and the particle 
position. The normalization of the results is analogous 
to the normalization used in Fig. 1 for ^pp. 



-pp _ 
Pi — 



12r;^P 



(45) 



except that the viscosity appears in the numerator, be- 
cause of the factor 77 in the definition (40) of /iPp. The 
normalized polarizability coefficient (45) is equivalent to 
the dimcnsionlcss dipole moment (44) produced by the 



pressure gradient ^p'^^'^ acting on a force- and torque-free 
particle. 

A comparison of the results shown in Fig. 2 with those 
depicted in Fig. 1 indicates that 

(46) 



Mr < cf^ 



where the equal sign holds for a particle touching the 
wall. In such configurations the particle cannot move, 
owing to the diverging lubrication forces. Thus the par- 
ticle polarizability is the same in the mobility and friction 
formulations. Relation (46) is consistent with a general 
observation that the energy dissipation is larger in a sys- 
tem with more constraints, which can be demonstrated 
using variational techniques [37]. Relation (46) can also 
be directly obtained from (40), by observing that the 
translation-rotation resistance matrix is positive definite. 

The results shown in Fig. 2 indicate that /i^'^ has a min- 
imum a,t Z = H/2, and that for large values oi H/{2a) the 
minimal value nearly vanishes. This is because for force- 
and torque-free particle and large wall separation, the 
dominant contribution to the dipole moment D comes 
from the stresslet induced on the particle due to the lo- 
cal velocity gradient of the imposed parabolic flow. At 
the center of the channel, the imposed velocity gradient 
vanishes, so there is no stresslet contribution. 

Our results for the translational and rotational polariz- 
ability coefficients /uP* and /xP"^ are shown in Figs. 3 and 4. 
To emphasize the effect of the particle size and position 
on the dipole moment D, the translational polarizability 
is normalized by the coefficient 



l^l\Z) 



w 



'Z{H-Z), 



(47) 



corresponding to the dipole moment (18) of a unit lateral 
point force (Stokeslet) applied to the fluid at the position 



Z . Similarly, the rotational polarizability is normalized 
by the coefBcient 



VU7\ - 



p'.iz) 



¥- 



\\H^Z\ 



(48) 



representing the dipole moment associated with the dipo- 
lar far-field flow produced by a unit lateral point torque 
(rotlet) [23, 25, 38]. Note that (47) is symmetric and (48) 
is antisymmetric with respect to the channel center. 

The results in Figs. 3 and 4 indicate that /iP*//ig ~ 
[i^^ I [i^ ~ 1 for H/{2a) ^ 1, except for the regions ad- 
jacent to the walls, where the normalized polarizabilities 
decrease logarithmically to zero for a particle touching 
a wall. The normalized polarizability /iP'/^o decreases 
monotonically with decreasing wall separation, whereas 
fjP'^ / jIq'^ is non-monotonic in H/{2a). 



V. MACROSCOPIC SUSPENSION FLOW 

In this section we apply the results of the above anal- 
ysis to determine the average volume flux in a dilute 
suspension bounded by two parallel planar walls. The 
average flow can be driven by the macroscopic pressure 
gradient, by lateral forces or by lateral torques applied 
to the particles. 

To determine macroscopic equations governing suspen- 
sion flow, we first show that the far-field Hele-Shaw fiow 
(1) driven by the dipolar pressure distribution (3) con- 
tributes to the macroscopic suspension velocity. This be- 
havior is analogous to the 2D electrostatics, where the 
electric field produced by induced dipoles contributes to 
the macroscopic electrostatic displacement field (as dis- 
cussed in Sec. II C 1). 

The electrostatic analogy can be directly applied to 
fluid flow through a system of immobile particles in a 
channel. In this case the fluid transport is governed by 
the Darcy's equation 



-K,V 



\\P, 



(49) 



which is a counterpart of the constitutive relation be- 
tween the electric and electrostatic-displacement fields. 
In equation (49) V||p denotes the macroscopic pressure 
gradient, and u is the macroscopic velocity, defined as 
fluid-flux density averaged across the channel. The effec- 
tive permeability coefficient k plays a role similar to that 
of the effective dielectric constant in the corresponding 
electrostatic system. Equation (49) can be obtained by 
combining Eq. (9) with the "pp" component of Eq. (28). 

Note that relations (9) and (49) involve the macro- 
scopic pressure gradient (rather than the gradient of the 
external pressure p™'), to ensure that the constitutive re- 
lations involve only local quantities and are independent 
of the boundary conditions. 

For a system of particles suspended in a fiuid, sus- 
pension transport is described by a constitutive equation 
that includes additional force and torque contributions. 



where fPP, i/P*, and i/P'' are the effective mobility coeffi- 
cients. 

To derive Eq. (9) and determine the transport coef- 
ficients in Eqs. (49) and (50), we consider the average 
pressure gradient and average velocity in a periodic sys- 
tem representing a macroscopically uniform quasi-two- 
dimensional medium. 



A. Periodic Green's function 

We start our analysis by deriving appropriate expres- 
sions for the periodic Green's functions for Stokes flow in 
a parallel-wall channel. In the present paper, our formu- 
lation is used as a theoretical tool to obtain the relation 
between the quantities D and u. Our explicit results are 
also applied in numerical calculations presented in Sec. 
VI. More generally, our formulas for the periodic Green's 
functions can be employed in Stokesian-dynamics and 
boundary-integral algorithms for parallel-wall geometry. 



1. Near-field and far-field contributions 

In our approach, the periodic flow and pressure Green's 
functions TP°'' and Qp°'' are evaluated by splitting them 
into the asymptotic Hele-Shaw parts T^g and Qh's ? ^^d 
quickly convergent lattice sums of 3D corrections to the 
asymptotic Hele-Shaw behavior, 



^(r,r') 



.'\ _ rppcr 
'-HS 



(r,r') + ^<5T(r,r:,), 



(51a) 



QP-(r, r') = QP's'-(r, r') + ^ SQ{r, r'J. (51b) 



Here n = {ux, %) (with nx,ny — 0, ±1 . . . ) are the in- 
dices of the periodic lattice. 



r^ = r' + rixL^ex + UyLye 



(52) 



TJu = lyp' • :r -f lyp^ . T - i^pp V||p, 



(50) 



are the positions of the periodic images of the source 
point r', and L^ and Ly are the lattice constants. The 
near-fleld contributions ST and (5Q are deflned as the dif- 
ferences between the exact and asymptotic non-periodic 
Green's functions, 

,5T(r, r') = T(r, r') - THs(r, r'), (53a) 



(5Q(r,r')-Q(r,r')-QHs(r,r'). (53b) 

Since the flow and pressure flelds in a parallel- wall chan- 
nel tend to the asymptotic Hele-Shaw form exponentially 
on the lengthscale H [23, 24], the near- field contributions 
(53) exponentially vanish at large lateral distances 5p be- 
tween the field and source points r and r'. Our numerical 
tests [23] indicate that the exact and asymptotic Green's 
functions are nearly identical for 5p/ H > 3. Therefore, 
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in practical calculations only a small number of terms 
need to be evaluated to determine the lattice sums in 
Eqs. (51) with high accuracy. 

We note that explicit expressions for all terms in Eqs. 
(53) are known [cf. Eq. (15) and the results in Appendix 
A]. Analytic formulas for the asymptotic Hele-Shaw 
Green's functions 



^HS ^^^ derived in the fol- 
lowing section and in Appendix C. 



TP-g-- and Q^ 



2. Far-field Green's functions for periodic system 



HS ^^"i Qhs 



To evaluate the far-field components T^g"^ 
of the periodic Green's functions (51) we start from the 
direct lattice sums 



TP-(r,r') = ETHs(r,r'„), 



-HS' 



Q&s(r,r') 



J2Qm{r,r'J. 



(54a) 



(54b) 



Combining relations (54) with (15) allows us to express 
the hydrodynamic periodic Green's functions in the Hele- 
Shaw regime in terms of the periodic solution of the cor- 
responding electrostatic problem. Accordingly, we have 



Tlis(r,r') 



iry-iz(/7-z)VQP^s^(r,r' 



(55a) 



QP's''(r, r') = -3TT-^H-^Vwip - p')z'{H - z'), (55b) 

where the Wigner potential w{p — p'), is the periodic 
solution of the Poisson equation 



\/'^w{p- p') = -27r 



E^(p 



(56) 



Here p^ is the lateral component of the lattice vector 
(52), and Sp = L^Ly is the area of a unit cell [39]. 

The Wigner potential w{p — p') can be determined 
using standard Ewald summation techniques [40, 41]. 
Well-developed accelerated algorithms for calculating 
this function are also available [40]. Equations (51) and 
(54) thus reduce the problem of evaluating the doubly- 
periodic 3D hydrodynamic Green's functions to a much 
simpler scalar problem (56). Explicit expressions for the 
Wigner function and its multipolar projections are given 
in Appendix C. We also derive there Ewald sums for the 
hydrodynamic Green's functions (55). 



B. Average flow produced by a point force 

The asymptotic periodic pressure Green's function 
(55b) is proportional to the gradient of the periodic 
Wigner function w. Thus Qh*s i^ normalized to yield 



zero average value over a unit cell (but any constant can 
be added without changing the system dynamics). 

In contrast, the velocity Green's function TP°'' is fully 
determined by Eqs. (51a) and (55), with no gauge con- 
stants involved. Adding a constant would violate the 
boundary conditions on the channel walls — thus, the av- 
erage flow field in a wall-bounded periodic system cannot 
be set independently. Instead, the average flow is a func- 
tion of the applied pressure drop and the induced-force 
distribution (unlike the corresponding behavior in the in- 
finite space). 

As a key step in our analysis of suspension flow in a 
channel, we derive an expression for the average flow Ug 
produced by a lateral point force F|| applied at a point 
r' = (p', z'). (Only a lateral force needs to be considered, 
because a transverse force F^ = F^e^ does not produce 
an average flow, by symmetry). The overall pressure drop 
in the system is assumed to vanish. 

The average flow is given by the integral over a unit 
cell 



Uo = f^p 



vo(r)dr 



of the flow field 



vo(r)==TP-(r,r')-F 



(57) 



(58) 



produced by the applied force. To determine the volume 
integral in Eq. (57) we use the decomposition (51a) of 
the Green's function TP°'^ into the near-field and far-field 
components. The far-field contribution vanishes because 



T&s(r,r') 



0, 



(59) 



owing to Eq. (55a) and the periodicity of the pressure 
Green's function Qh"s • Therefore, the average flow Uq is 
associated with the integral of the near-field contribution 



<5vo(r)=^<5T(r,r'„ 



II 



(60) 



over the unit cell Op. Using the invariance of (5T(r,r') 
with respect to lateral translations, this integral can be 
represented as 



Uo 



VLr 



(5vo(r) dr = f7 



p-i/ (5T(r,r')-F|| 



(61) 
where i7oo — Tlm^vi^) i^ the whole infinite volume of 
the channel. 

The integral (61) can be related to the dipolar strength 
(18) of the asymptotic pressure distribution (16) by in- 
serting under the integration sign the identity tensor 
I ~ Vr and integrating by parts. Since the boundary 
term vanishes, owing to the boundary conditions on the 
walls and the rapid decay of 5T(r, r') for p ^- oo, we find 
that 



Uo 



-VLp-^ I rV-(5T(r,r')-F| 



dr. 



(62) 
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Noting that the exact Green's function in Eq. (53a) is 
divergence-free and using relations (15), (16), and (20) 
for the asymptotic contribution we obtain the following 
fundamental result 



uo = Koflp ^iJDo, 



(63) 



where Dq is the dipole moment (18), and kq is the per- 
meability coefficient for the particle-free channel (10). In 
an explicit form we have 



Uo 



^T]'^z'{H -z')np-^F 



(64) 



Relation (64) is consistent with the average fluid velocity 
produced by a planar force distribution 



fsir) 



np-^ Hd{z 



z')Fu 



(65) 



C. Average flow in the particle presence 



The macroscopic fields V||p and u that appear in the 
effective-medium equations (49) and (50) can be identi- 
fied with the volume averages of the microscopic pressure 
gradient V||p and velocity v over a unit cell ftp, 



V||p = Up 



Qr, 



V||pdr, 



vdr. 



(66a) 



(66b) 



To determine the macroscopic transport coefficients in 
equations (49) and (50), we thus need to derive appro- 
priate expressions for these averages. 

In this section we consider a system of N particles at 
positions Hi (i ~ 1, . • . , N) in the unit cell. The parti- 
cles are represented by the corresponding induced-force 
distributions Fi(r). 

It is convenient to represent the pressure and flow fields 
(12) (generalized to a multiparticle system) as the super- 
positions of the external and scattered contributions 



N 



N 



Y.<- 



(67a) 



(67b) 



We assume that the applied pressure gradient Vp"''* is 
constant in space and has only lateral components x and 
y. The corresponding external velocity field (27) depends 
only on the transverse coordinate z. The scattered flow 
and pressure fields v' and p' are periodic, and they are 
given by the expressions 



V^ 



QP^'(r,r')-F,(r')dr', 



(68a) 



^(r,r')-F,(r')dr'. 



(68b) 



By integrating the lateral gradient of the pressure (67a) 



over the unit cell fJp we find 



V||p = Vf 



(69) 



where the integrals of the perturbation-pressure terms 
V||p^ vanish by periodicity of p'. Equation (69) is im- 
portant, because it allows us to express the macroscopic 
flow and particle motion in a channel in terms of the 
macroscopic pressure (rather than the external pressure) 
in the constitutive relations derived in Sec. VI. 

Integrating relation (67b) with the external flow given 
by Eq. (27) yields 



u = -Ko Vp^"* + rip 



N 



v^(r)dr. (70) 



The integral on the right-hand-side of the above equation 
can be evaluated using the result (63) for the average ve- 
locity produced by a point force. Combining point-force 
results (57) and (58) with (68b), and taking into account 
that the transverse force components do not contribute 
to the average flow, we find 



N 



U = Koi-Vp'^^' +np-^HJ2^^), 



(71) 



i=l 



where D^ is the dipole moment (22) of particle i. Defining 
the average dipole moment 



N 



D = iV^^D„ 



(72) 



relation (9) for the average flow is thus obtained. 

According to the above derivation, equation (9) is valid 
for arbitrary particle densities, provided that the induced 
forces Fi are evaluated with the multiparticle hydrody- 
namic interactions properly taken into account. In the 
low-density limit, the average dipole moment (72) can be 
expressed in terms of the dipole moment of an isolated 
particle D(Z), 



risD 



H-a 



n{Z)T>{Z)AZ, 



(73) 



where n{Z) is the local particle number density per unit 
volume, averaged over the lateral position within a unit 
cell. 



VI. TRANSPORT COEFFICIENTS 

In this section we use the relation for the average ve- 
locity (9) to obtain the effective macroscopic equations 
for fluid and particle transport in a parallel-wall channel. 
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FIG. 5: Effective permeability (76) of a channel with a reg- 
ular particle array adsorbed on the lower wall, versus area 
fraction cj)s, for several values of normalized channel width 
H /{2a) (as labeled). Hexagonal particle arrangement (dotted 
lines); square arrangement (dashed lines). 



FIG. 6: Effective permeability (76) of a channel with a hexag- 
onal particle array adsorbed on the lower wall, versus area 
fraction (J)b , for two values of normalized channel width H /{2a) 
(as labeled). Exact result (solid lines); low-density limit (78) 
(dotted); Clausius-Mossotti approximation (79) (dashed). 



In a dilute-suspension regime, the macroscopic transport 
equations are obtained by combining (9) with the fric- 
tion relation (28) or mobihty relation (36). At higher 
particle concentrations, we use (9) in combination with 
linear constitutive relations between the dipole moment 
and the macroscopic forcing for a system of interacting 
particles. 

We consider here two important cases. In Sec. VI A we 
discuss fluid transport through a fixed bed of particles, 
and in Sec. VI B we examine transport of a suspension of 
freely moving particles. 



A. Fixed bed of particles 

1. Permeability coefficient 

We begin our analysis of fluid transport through a fixed 
particle array by considering a low-density system. In 
the low-density domain, the dipole moment D contribu- 
tion to the average flow (9) can be obtained by averaging 
the polarizability component of the generalized resistance 
relation (28) over the particle distribution. Taking into 
account that U = Jl = we find that 



^ff^D 



•77 



IIP, 



where 



*pp 



H-a 



n{Z)CP{Z)dZ 



(74) 



(75) 



is the mean value of the ^pp component of the gener- 
alized resistance matrix in Eq. (28). The driving force 
in the polarizability relation (74) is the gradient of the 
macroscopic pressure Viip, which at low densities is iden- 
tical to V||p°''*. By inserting (74) into (9) we obtain 



Darcy's equation (49) with the permeability coefficient 
of the form 



K„{1 - 12n,H-^rj-^en, 



(76) 



where definition (10) was used to factor out the perme- 
ability of particle- free channel kq . 

The polarizability relation (74) and expression (76) 
for the channel permeability are valid not only for 
dilute-suspensions, but also for arbitrary particle con- 
centrations, provided that the average (75) is replaced 
by the corresponding relation applicable in the high- 
concentration regime. As in other problems of field prop- 
agation through random media, the polarizability coeffi- 
cient (^PP can be expressed in terms of cluster integrals 
that involve pair, triplet, and higher-order resistance 
functions. Derivation of such a relation requires apply- 
ing an appropriate renormalization procedure in which 
the external pressure is replaced by the macroscopic pres- 
sure in order to obtain absolutely convergent results for 
^PP. Alternatively, we can get the polarizability from the 
dipole moment of particles in a channel with periodic 
boundary conditions in the lateral directions. Numeri- 
cal results described in the following section have been 
obtained using the latter technique. 



2. Numerical results 

To illustrate the effect of immobile particles on the per- 
meability of a parallel- wall channel we present results for 
a particle monolayer adsorbed on one of the walls. Fig- 
ures 5-8 show the permeability coefficient (normalized by 
the permeability of a particle-free channel) for hexagonal 
and square particle arrays. For the square and hexag- 
onal symmetry, the tensorial permeability coefficient is 
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FIG. 7; Effective permeability (76) of a channel with a 
closed-packed regular particle array adsorbed on the lower 
wall, versus normalized channel width H/{2a). Hexago- 
nal particle arrangement (dotted lines); square arrangement 
(dashed lines). Inset shows a blowup of the plot for small 
values of H /(2a). 



isotropic 



K = k\\\ 



(77) 



where l|| is the lateral unit tensor (21). 

Figure 5 illustrates the dependence of the permeability 



coefficient k on the particle area fraction 



for 



several values of channel width. The results indicate that 
the permeabilities of the hexagonal and square arrays are 
nearly identical functions of 4>s — only for a channel with 
H/{2a) K, 1 there is a noticeable difference, especially at 
high surface coverage 0s- 

In Fig. 6 our numerical results for hexagonal arrays are 
compared to the low-density limiting behavior 



- = i-Cf 



KO 



(78) 



where the normalized one-particle polarizability coeffi- 
cient Cf^ is given by Eq. (43). Relation (78) follows from 
Eqs. (75) and (76) applied to a particle monolayer. For 
particles positioned in the midplane of the channel the 
particle density, averaged over the lateral directions, is 
n{Z) = ns5{Z — \H), where fig = N/{LxLy), and N de- 
notes the number of particles in a unit cell. For more 
general particle distributions, the average value of (^f^ 
would appear in (78). 

In Fig. 6 we also plot the Clausius-Mossotti approxi- 
mation 



1 - irp 



PP. 



Ko 



1 + Ki^ 



pp. 



(79) 



which is a generalization of the classical electrostatic 
Clausius-Mossotti formula [28] to our present problem. 
The results indicate that for tightly confined systems 



FIG. 8: Fluid-flux ratio (82) for the same two systems as 
those represented in Fig. 7. 



with H/{2a) ~ 1, the approximation (79) is quite accu- 
rate in the area-fraction range 4>s ^ 0.2. For weaker con- 
finements, the range of validity of the Clausius-Mossotti 
formula is smaller. 

The permeability of close-packed arrays (with the 
close-packing area fraction (j>s = \/3tt/Q for hexagonal 
and 03 — 7r/4 for square ordering) is plotted in Fig. 7 
versus the dimensionless wall separation H/{2a). The 
results indicate that for H/{2a) w 1, the permeability 
coefficient is reduced to about 1 % of the permeability 
kq of a particle-free channel with the same width. Such 
a significant reduction of the channel permeability was 
observed in recent experiments [42]. 

At weaker confinements, the reduction of the perme- 
ability is much smaller, but it is still quite significant for 
H/{2a) < 10. For moderate and large channel widths 
the hindrance of fluid flow can be accurately accounted 
for in terms of the reduced effective channel width. In 
this approximation, the particle array is replaced by an 
equivalent solid slab of width Zns , occupying the region 
< z < Z^s- The permeability of the narrowed channel 



where 



Kcff 



H, 



T2^ 



H 



-1tt2 



^NS 



(80) 



(81) 



is the effective channel width. 

The accuracy of the effective-width approximation (80) 
can be estimated from a plot of the ratio 



A 



Hk 



HcS Kcff 



(82) 



of the fluid flux through the channel with adsorbed parti- 
cles to the fluid flux through a particle- free channel of the 
reduced width. Such a plot is shown in Fig. 8 for closely 
packed hexagonal and square arrays. From our data we 
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find Zns = 0.907 for the hexagonal particle array and 
■^NS = 0.895 for the square array. With the above values 
we get A w 1 for H/{2a) > 2. In this range of channel 
widths, the approximation that neglects the roughness 
of a densely packed particle layer adsorbed on the wall 
is thus accurate. We note that our findings are consis- 
tent with earlier investigations of fluid flow near rough 
surfaces [43, 44]. 



B. Suspension of freely moving particles 

1. Macroscopic constitutive equation 

The linear constitutive equation relating the macro- 
scopic fluxes to macroscopic forces in the dflute- 
suspension regime is obtained by combining the dipolar 
expression for the average suspension velocity (9) with 
the generalized mobility relation (36). Relation (36) is 
first averaged over the particle distribution, which yields 
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(83) 



(84) 



similar to Eq. (75) in the friction-representation case. In 
Eq. (83) it is assumed that the same force ^ and torque 
T act on all the particles in the system. We also assume 
that both J^ and T have only the lateral components. As 
it has been done for immobile particles, in Eq. (83) the 
external-pressure gradient is replaced with the gradient 
of the macroscopic pressure Viip to obtain a local consti- 
tutive relation that can be generalized to dense systems. 
By combining (9) and (83) we find 
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In the constitutive relation (85), the macroscopic fiuxes 
are: the particle flux fisU, particle angular flux ngU, and 
the suspension volume flux Hu. 

Similar to the immobile-particle case, relations (83), 
(85), and (86) are valid at arbitrary densities, provided 
that the low-density formula (84) for the generalized mo- 
bility coefficients p, is replaced with the corresponding 
expression that is appropriate at high densities. 
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FIG. 9: Effective permeability (86b) of a channel with a reg- 
ular particle array suspended in the midplane Z — H/2, ver- 
sus area fraction (jf>s, for several values of normalized channel 
width H/{2a) (as labeled). Hexagonal particle arrangement 
(dotted lines); square arrangement (dashed lines). 



2. Onsager reciprocal relations for the generalized mobility 
matrix v 

In Sec. IV we have shown that at low particle concen- 
trations the generalized mobility/polarizability matrix /x 
in Eq. (36) satisfies the Lorentz symmetry (8). The same 
symmetry also applies to the average defined by Eq. (84). 
It follows that the matrix of the kinetic coefficients (86) 
in the constitutive relation (85) satisfies the Onsager re- 
ciprocal relation 



.AB 



,SAt 



(87) 



While our derivation is given here only for dilute sus- 
pensions, one can show that the symmetry relation (87) 
is valid at arbitrary concentrations. The symmetry of the 
matrix u can be demonstrated using arguments similar 
to the ones given in Sec. IV, but applied to a periodic 
multiparticle system. 



3. Numerical results 

At low suspension concentrations and, more generally, 
for suspensions that are isotropic in the lateral direc- 
tions, the transport coefficients u^^ are proportional to 
2D isotropic tensors 
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and 



i^^' = v^ 



(88b) 



where e|| is the lateral alternating tensor (42). Relations 
(88) are also satisfied for hexagonal and square particle 
lattices. 
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FIG. 10: Effective permeability (86b) of a channel with a 
hexagonal particle array freely suspended in the midplane 
Z — H/2, versus area fraction (j)s, for two values of normal- 
ized channel width H/{2a) (as labeled). Exact result (solid 
lines); low-density limit (89) (dotted); Clausius-Mossotti ap- 
proximation (90) (dashed). 



The effective channel permeability coefficient j^pp, nor- 
malized by the permeability of a particle-free channel, is 
plotted in Fig. 9 for infinite hexagonal and square parti- 
cle arrays moving in the midplane of the channel. As for 
arrays of immobile particles, we find that the permeabil- 
ities of the square and hexagonal arrays with the same 
area fraction are nearly the same. 

A comparison of the results depicted in Figs. 5 and 9 
indicates that the permeability of a channel with parti- 
cles freely suspended in the midplane Z = H/2 is much 
higher than the corresponding permeability of a channel 
where the particles are adsorbed at a wall. According to 
Fig. 9 the correction to the permeability due to the parti- 
cle presence is below 10% in the regime H/{2a) > 2, even 
for close-packed arrays. This result is consistent with the 
small value of the single-particle polarizability for a par- 
ticle at the symmetry plane Z = H/2, as illustrated in 
Fig. 2. Both the small particle polarizability and the 
minor correction to the channel permeability stem from 
the small velocity gradient in the midplane of the chan- 
nel, which implies that the suspended particles do not 
significantly perturb the fluid flow. 

In Fig. 10 the results of the numerical calculations 
for periodic particle arrays are compared with the low- 
density expansion 
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where jl^ is given by (45). We also plot the Clausius- 
Mossotti approximation 
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FIG. 11; Effective mobility coefficient v"^^ for an infinite par- 
ticle array in the midplane Z — H/2, normalized by the corre- 
sponding result for an array of point forces (91), versus area 



fraction 



for several values of normalized channel width 



H/(2a) (as labeled). Hexagonal particle arrangement (dotted 
lines); square arrangement (dashed lines). The coefficient v'^^ 
describes the average suspension flow produced by a lateral 
force applied to the particle array. 



dicatc that the Clausius-Mossotti approximation has a 
broader range of validity for freely suspended particles 
than for the adsorbed ones. 

The transport coefficients i^P* and i^p'', representing the 
average flow produced by a lateral force and torque act- 
ing on the particles, are plotted in Figs. 11 and 12. The 
transport coefficients are shown normalized by the re- 
spective results 



pt/ 
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A comparison of the results shown in Figs. 6 and 10 in- 



for the average flow produced by arrays of point forces 
and torques applied to the suspending fluid at the particle 
positions Z, where /Iq and /j.q'^ are given by relations (47) 
and (48). The results are presented for square and hexag- 
onal particle arrays. For the system driven by the lateral 
force, the particles are in the midplane of the channel 
Z = H/2. For the lateral torque the rotational mobility 
coefficient i^P'' is shown for Z = a -I- ^{H — 2a), because 
the torque applied to the particles produces average flow 
only for off-center positions. 

The results in Fig. 11 indicate that for H /{2a) > 2 the 
average flow produced in a channel by an external force 
applied to the particles is well represented by the point- 
force approximation (91), even for dense particle arrays. 
The normalized average flow in this regime is insensitive 
both to the dimensionless channel width H/{2a) and to 
the area fraction <j>s. For smaller values of H/{2a), the 
normalized average flow is smaller than the flow produced 
by point forces. However, a significant deviation of the 
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FIG. 12; Effective mobility coefficient v^'^ for an infinite par- 
ticle array in the plane Z — a + ^{H — 2a), normalized by the 
corresponding results for an array of point torques (92), ver- 
sus area fraction (jia, for several values of normalized channel 
width H/{2a) (as labeled). Hexagonal particle arrangement 
(dotted lines); square arrangement (dashed lines). The coef- 
ficient u^" describes the average suspension ffow produced by 
a lateral torque applied to the particle array. 



normalized transport coefficient v^^ jv^ from unity is ob- 
served only for particles with the diameter nearly equal 
to the channel width. This behavior is similar to the one 
seen in Fig. 9 for arrays driven by an external pressure 
gradient. The low-density values are consistent with the 
one-particle results represented in Fig. 3. The results for 
hexagonal and square lattices for a given area fraction 
are nearly indistinguishable. 

As depicted in Fig. 12, the average flow produced by 
a torque T applied to the particles strongly depends on 
the particle area fraction, even for large wall separations. 
This is because particle rotation involves relative mo- 
tion of the surfaces of the spheres, whereas the motion 
of force-driven particle monolayers does not involve any 
relative particle displacements. Moreover, we find that 
the results for the square and hexagonal particle arrays 
are significantly different for (/)s ^ 0.5. The transport 
coefficient v^'^ vanishes at the close-packing area frac- 
tion for a given system geometry, because the particle 
rotation is arrested by the lubrication forces between the 
touching particles. In the low-density regime the results 
are independent of the particle lattice. The low-density 
values of the transport coefficients are equivalent to the 
one-particle results represented in Fig. 4. 



VII. CONCLUSIONS 

We have presented a detailed analysis of the far-field 
scattered flow produced by spherical particles in Stokes 
flow bounded by two parallel planar walls. Furthermore, 
we have examined the effect of the far-fleld particle re- 



sponse to external forcing on the macroscopic suspension 
flow. Both the permeability of a system of flxed particles 
in a channel and the macroscopic dynamics of a suspen- 
sion of freely moving particles were analyzed. (We note 
that related ideas were also explored in [45], in the con- 
text of molecular-dynamics simulations of the motion of 
nano-particles in a confined fiuid.) 

For a system of fixed particles, the macroscopic fluid 
flux is related to the macroscopic pressure via linear 
Darcy's law. We have shown that for a given macroscopic 
pressure gradient, the difference between the fluid flux in 
the particle presence and in a particle-free channel can 
be expressed in terms of the effective 2D dipole moment 
characterizing the amplitude of the far-field Hele-Shaw 
dipolar scattered fiow produced by the particles. From 
this amplitude we have evaluated the particle contribu- 
tion to the effective permeability coefficient. 

A similar physical picture also applies to the fiow of 
a suspension in a channel. However, in this case the 
macroscopic volume flux can be produced not only by 
the macroscopic pressure gradient but also by the ex- 
ternal force or torque applied to the particles. There 
are also particle fluxes corresponding to the linear and 
angular particle velocities. Therefore, the constitutive 
relation for suspension fiow through a channel involves 
three forcing and three flux components. The macro- 
scopic fluxes and driving force are related through a 3 x 3 
matrix of (generally tensorial) transport coefficients. We 
have demonstrated that, with a proper normalization, 
this matrix is symmetric, i.e., the transport coefficients 
satisfy the Onsager reciprocal relations. 

Our theoretical analysis has been supplemented by nu- 
merical results for transport coefficients describing dy- 
namics of square and hexagonal particle arrays (for par- 
ticle monolayers adsorbed on a wall and monolayers of 
freely suspended particles). We have shown that dense 
arrays of tightly confined particles can reduce fluid flow 
through a channel by as much as 99%. 

We have also proposed generalized Clausius-Mossotti 
formulas for the channel permeability, both for fixed and 
freely suspended particles. These formulas are analo- 
gous to the well-known electrostatic Clausius-Mossotti 
approximation. At moderate particle concentrations, our 
expressions agree well with the numerical calculations, 
especially for tightly confined particle arrays. 

Our numerical results for the transport coefficients 
characterizing macroscopic motion of regular particle ar- 
rays can be used to describe the macroscopic deforma- 
tion of finite-size regular 2D particle clusters. As we have 
shown in our recent paper [20] , such arrays exhibit a com- 
plex nonlinear dynamics that involves rearrangements of 
a deformed particle lattice. There are also order-disorder 
transitions resulting from lattice instabilities. In future 
publications we will analyze these problems using the 
macroscopic theory developed in this paper. We will also 
determine particle and fluid transport in suspensions of 
randomly distributed particles. 

Our analysis of suspension flow in parallel-wall chan- 
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nels with periodic boundary conditions has also pro- 
vided another important result: wc have derived explicit 
Ewald-summation formulas for the flow and pressure pe- 
riodic Green's functions in the parallel-wall geometry. 
These formulas can be applied in Stokesian-dynamics and 
boundary-integral simulations of suspension and emul- 
sion flows in narrow channels and slit pores, and in stud- 
ies of the dynamics of confined macromolecules. 

In future publications we will discuss the effect of the 
far-field flow on the dynamics of macromolecules (e.g., 
polymer chains). Earlier studies have suggested that in 
a confined system the far-field contribution to the hy- 
drodynamic interactions between the chain segments can 
be neglected, because, on average, this contribution van- 
ishes by symmetry [46, 47]. Our preliminary investiga- 
tion based on the present results indicates that the role of 
the far-field fiow is much more subtle, although the far- 
field hydrodynamic interactions do not change the Rouse 
scaling exponent for the longest relaxation time (only the 
prcfactor is affected). 
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pressions are obtained using our Cartesian-representation 
approach introduced in [3, 4]. In our previous papers ex- 
plicit results were given only for multipolar projections 
of the Green's tensor T. 

The Green's functions for Stokes flow between two par- 
allel walls can be expressed as a sum of the free-space part 
and the wall contribution. 



T(r,r') = To(r-r')+T'(r,r') 



Q(r,r') = Qo(r-r') + Q'(r,r') 



where 



To(r) = -^(i + rr), Qo(r) = 
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(Ala) 



(Alb) 



(A2) 



are the Oseen tensor and the corresponding pressure 
Green's function. As in [3, 4] the wall contributions to the 
Green's functions, T' and Q', are represented in terms 
of lateral Fourier integrals of simple matrix products. 

It is convenient to express components of the tensor 
T' and vector Q' in terms of the spherical basis of unit 
vectors [481 



1 



e-i 



(err - iey), eo 



V2' 
Accordingly, we have 



ei 



1 

'V2 



\^x I IG^ 



T'(ri,r2)= Y^ Yl rmm'(i"i,r2)e,„e; 



(A3) 



(A4a) 



m— — l rn,' — — 1 



APPENDIX A: VELOCITY AND PRESSURE 
GREEN'S FUNCTIONS T AND Q 

In this Appendix we present our explicit expressions 
for the velocity and pressure Green's functions T and Q 
for Stokes flow between two parallel planar walls. The ex- 
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The components T^^, and Q^, can be evaluated from 
the following 2D Fourier integrals 
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(A5b) 



where pi2 = Pi — P2, and k — (k, tp) represents the wave 
vector k in polar coordinates. The integral kernels in 
Eqs. (A5) can be expressed as products of several simple 



matrices. 



tmm'ik;zi,Z2) = B''' (zi , fc, m) • Ztw (k) • B(z2,fc,m'), 

(A6a) 

qm'ik;zi,Z2) = b'''(zi,fc)-ZTw(k) • B(z2,/c,m'), (A6b) 
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where the dagger denotes the transpose. 
As explained in [3, 4], the matrix 



ZTw(k) 



Z-i S++{-kH) 
Sc-{kH) Z-i 



(A7) 



describes the multiple reflections of Cartesian hydrody- 
namic basis fields from the parallel walls. The 3x3 
component displacement matrices 



S++{-kH) = [S^-(kH)]^ 
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describe the propagation of the flow flelds between the 
walls, and the matrices 
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represent scattering of the flow field from the walls. The 
matrices 



B{k,z,m) = [2(l-m)!(l + m 
where 


),]-l/2 


■ B^{k,z,m) 
BU(fc,z,m) 


VIO) 


B^{k,z,m)^ (-l)"+ic-'=^ 


" -2kz + 2m^-l' 

2m 
1 


; 




± 

(Alia) 


BU(fc,z,m) = e-'=(^-^) 

describe the expansion of 
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correspond to the pressure (at point z), associated with 
the Cartesian basis fields centered at the lower or upper 
wall. The three components in the matrices (AS), (A9), 



(All), and (A13) correspond to the pressure, vorticity, 
and potential basis solutions of Stokes equations. 

The 2D Fourier integrals (A5) can be converted into 
the ID Hankel transforms by performing the angular in- 
tegration with the help of the relation 

(.27r 

(27r)-i / e''"'^e"'"^iMV' = i'"e''"'^i=J™(fcpi2), (A14) 



where (pi2 is the polar angle of the vector pi2, and 3m{x) 
is the Bessel function of the order m. The resulting ex- 
pressions are 
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(A16b) 
Substituting (A15) into (A4) and using definitions 
(A3) of the spherical basis vectors we get 

T'(ri,r2) = tii(ri,r2)l||+ti_i(ri,r2)(l|| -2P12P12) 

+V2[ioi (ri , r2)e2pi2 + ho (ri , r2)pi2ez] 
+too(i"i,i"2)e2e2, (Al7a) 

Q^,(ri,r2) = \/2(7i(ri,r2)pi2 + <Zo(ri,r2)ez. (A17b) 

We note that relations (A4)-(A17) for the Green's 
functions T and Q are equivalent to those derived by 
Jones [2], but our expressions are more transparent. In 
particular, the dependence of the integrands on the vari- 
ables zi and Z2 is clearly factored out, because z appears 
only in the matrices (AlO) and (A12). Moreover, our ex- 
pressions can easily be adapted to other boundary condi- 
tions at the walls (e.g., a fluid-fluid interface) by simply 
replacing the reflection matrix Z^. 



APPENDIX B: TRANSFORMATION VECTORS 

X(^ I Inia) AND X{lma \ A) 

In this Appendix we give explicit expressions for the 
transformational vectors X(A | Ima) (where A = t,r,p), 
deflned in Eqs. (30) and (31). Formulas for transla- 
tional and rotational transformation vectors X(t | Ima), 
'K.ilma I t), X(r | Ima), and 'K.{lma \ r) were derived in 
[4]. The transformation vector yi(lma \ p) (in a slightly 
different notation) is given in [5]. The remaining trans- 
formation vector X(p | Ima) can be directly obtained 
from relations (50)-(52) in Ref. [23]. 
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The transformation vectors X in Eqs. (30) are given 
by the expressions 



X(t I Ima) = {Uf/HaS^^i 



(Bla) 



formation vectors '^^{Iracr \ A) are related to (Bl) via the 
symmetry (32). 

The transformation vector X(p | Ima) is nonzero only 
for m = ±1 and 



X(r I Ima) = -2^Ti)^'^5n5„ier. 



(Bib) 



X(p I Ima) ^ 2"i/VC(Z; lma)em, (Blc) 

where e^ denote the basis vectors (A3), and m = 
— 1,0, 1. For other values of m, the transformation vec- 
tors vanish. It can be verified that the reciprocal trans- 
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The coefficient C in relation (Blc) can be expressed in 
the form [23] 
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where B^^{1; Z) denote the elements of the 3x3 matrix 
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The range A = 0,1,2 of the index A = Z — 1 in equation 
(B4) results from the conditions I > 1 and (B2). 



1. Multipolar periodic solutions of 2D Laplace 
equation 



APPENDIX C: FAR-FIELD CONTRIBUTIONS 
TO PERIODIC GREEN'S FUNCTIONS 

In this Appendix we list explicit formulas for the 
far-field contributions to periodic Green's functions for 
Stokes flow between two parallel walls. In Sec. C 1 we first 
consider the 2D scalar problem. We provide an Ewald- 
sum representation for the Wigner function w and list 
the corresponding expressions for the periodic multipo- 
lar solutions of Laplace equation in 2D. In Sec. C2 we 
give explicit formulas for the periodic Green's function 
(55) and for their spherical matrix elements (25). 
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The relations presented in this sections arc based on 
the results of Cichocki and Felderhof [41] who have de- 
rived Ewald-type expressions for the 2D periodic multi- 
polar potentials. 



a. Wigner function 

As shown in [41], Wigner function can be represented 
by the following formula 
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Here 
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are the direct and reciprocal lattice vectors, 

roc ^-xt 

Ei{x)= —j-At 



(C3) 



is the exponential function, and a is the splitting pa- 
rameter that controls the convergence of the direct and 



reciprocal sums. The gauge constant C„, — 1.3105329259 
[41] is used to set the hmit 



hm [w{p) - $0 (P - Pn)] = 0. 



b. Multipolar solutions 
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By definition, for ???, = we simply have 
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The periodic multipolar solutions $^ of the 2D 
Laplace equation are defined in terms of the non-periodic 
multipolar basis fields 
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For nonzero values of tti, the multipolar solutions are 
given by the expression 
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where r(j, x) is the incomplete Gamma function. 

Near the singularities at the lattice points p = Pn the 
periodic functions (C7) behave as 



*m(p)-*r„(p-pn)- 



(C8) 



For |m| > the gauge constants are determined by the 
condition that the fields $^ (p) can be expressed as com- 
binations of derivatives of the Wigner function w{p). 



c. Displacement theorems 

The multipolar fields (C7) can be used to construct the 
displacement formula for the Wigner potential [41] 



w{p + p) 



ZLxl^y 



E KM^i:ip)^ (C9) 



where the first term on the right-hand side corresponds 
to the background term in the periodic Poisson equation 
(56). Setting p = p' pi2 and using relations 



<i>+:(p-p')= E (-i)"g(rr.mo ^'7',+ '";,7 <i>±^(p)<i>+r(po 



m.'m 
m+m'— mi 



|?7l|! \m'\ 



(CIO) 



r 



[where 0{x) is the Heaviside step function] and we find a symmetric displacement relation 

2p.p'^<i>+(p)a>+*(p') + $li(p)$lr(p') (cii) 
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w{p + p,,-p') = -^{p' + p'^)+ J2 E <J'+(p)^+r(pi2;™|m')'i>+*(p'), (C12) 



^ m— — OO m' — — OO 



where 



IT f {\m\ -\- l/77'IV ~ 

5+7(pi2;m I m') = ----_(<5i„5w + <5-i™<5-i™') + ^(-"^™')(-l)" ^^^^^^$-,_„,(pi2) (C13) 



I 

is the displacement matrix for periodic multipolar scalar 2. Asymptotic Green's functions for Stokes flow 

fields. Integrating the above expressions with the multi- 
polar source distribution of order to', centered at p' = 

yields the displacement theorem for the multipolar peri- °- Green s functions Tjjg and Qjjg 

odic solutions (C7), 



TT 



00 



The relations given in the previous section can be used 
Kn'iP+Pi) = ^mO „'\ Pi+ E ^m{Pi)ScyiiPi'fn I ™'): *° dcrivc exphcit expressions for the asymptotic Green's 
ILxLy ,„^_oo functions (55). Taking the gradients of relation (C12) 

(C14) with respect to variables p and p', and evaluating the 
where the convergence condition pi < p is assumed. results at p = p' = we find 



Ths('^1''^2) = jT^zi{H - zi)z2{H - Z2) ^ ^ mm'S'+ 1 (P12; TO I TO')eme*„, , (C15a) 

m= — 1,1 m' = — 1,1 

QP^s^(ri,r2) = -^Z2(if-Z2) ^ ™'^+r(pi2; | ™')e:„. • (C15b) 

rn— — 1,1 



r 



The sums in (C15) can be evaluated explicitly using ex- 
pressions (C7), (C13), and (15), 



rppcr / N «J 



(l-l, r2) = 2 H^L L ""^^^ ^ ^l)^2(i? - ^2) l|| + E ^(2, TTCr ^jp^^ _ p„nTHs(pl2 - Pn; 21, Z2) 



HS 



47r 



T T y]^n'^Hs(k„;2;i,2;2)exp(-7rCT2fc^ + 27ripi2-k„), (C16a) 



QPf (n, r2) = Y^ r(l, 7rcr-2|p^2-p„nQHs(pi2-Pn; 22) + —^ ^ QHs(k„; ^2) exp(-7r(T2fc2+27ripi2 • k„), (C16b) 



" n/0 



where we have introduced notation The quickly convergent Ewald sums (C16) can be used for 

„ , , „ , , , „ s efficient evaluation of the periodic Green's functions T™ 

THs(pi2;2i,Z2) = THs(ri,r2), (C17a) — ^ "S 

applications. 

Qhs(Pi2; 21,22) = QHs(ri,r2)- (C17b) 



efficient evaluation of the periodic Green's functions T^g"^ 
and Qy" in Stokesian-dynamics and boundary-integral 



22 



Matrix elements 



The projections 
G^^ilma-.Yi I I'm'a' ■,r2) 



(w,+„.(ri) I TP^s I w 



I'm'a 



-(r2)), 
(C18) 



of the periodic Hele-Shaw Green's function T^g"^ onto the 
3D spherical basis can be obtained using relation (C12) 
and applying the method described in [49]. The results 
can be written in the form analogous to Eq. (46) in Ref. 

[23], 



G^ailma; ri 



I'm'a' ■,r2) 



irrjH^ 



C{Zi;lma)S+-{{pi2]m \ m')C{Z2]l'm'a'). 



(C19) 
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